Vacuum-induced coherence in quantum dot systems 
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We present a theoretical study of vacuum-induced coherence in a pair of vertically stacked semi- 
conductor quantum dots. The process consists in a coherent excitation transfer from a single-exciton 
state localized in one dot to a delocalized state in which the exciton occupation gets trapped. We 
study the influence of the factors characteristic of quantum dot systems (as opposed to natural 
atoms): energy mismatch, coupling between the single exciton states localized in different dots, 
different and non-parallel dipoles due to subband mixing, as well as coupling to phonons. We show 
that the destructive effect of the energy mismatch can be overcome by an appropriate interplay of 
the dipole moments and coupling between the dots which allows one to observe the trapping effect 
even in a structure with technologically realistic energy splitting on the order of milli-electron- Volts. 
We also analyze the impact of phonon dynamics on the occupation trapping and show that phonon 
effects are suppressed in a certain range of system parameters. This analysis shows that the vacuum 
induced coherence effect and the associated long-living trapped excitonic population can be achieved 
in quantum dots. 

PACS numbers: 78.67.Hc,78.47.Cd,03.65.Yz 
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I. INTRODUCTION 

Pairs of closely stacked quantum dots (QDs) coupled 
via inter-band dipole moments^ (double quantum dots, 
DQDs) or by tunneling resulting from carrier wave func- 
tion overlap and Coulomb correlations^— (quantum dot 
molecules, QDMs) attract much attention due to the 
richness of their physical properties which show huge 
technological promise for nanoelectronics, spintronics 
and quantum information processing applications. The 
unique features of these systems, as compared to in- 
dividual QDs, can be used as the basis for long-time 
storage of quantum information! 7 - and conditional opti- 
cal control of carrier states 8 - which pave the way to an 
implementation of a two-qubit quantum gateS. Double 
dot structures enable also coherent optical spin control 
and entangling 10 - 12 or may act as sources of entangled 
photons^. It has been shown that the exciton spectrum 
of a QDM can be used to define an excitonic qubit with 
an extended life time^, that information can be written 
on the spin state of the dopant Mn ion located in one of 
the dots forming a DQDA 5 -, and that a photon emitted 
by a nearby quantum point contact may induce carrier 
transfer in a DQD system^. The richness and complex- 
ity of the physical properties of these systems have been 
manifested in many optical experiments^ 7 - - — . 

The spacing between the dots forming intentionally 
manufactured DQDs and QDMs is typically on the order 
of nanometers which is two to three orders of magnitude 
smaller than the wavelength of a resonantly coupled pho- 
ton. Under such conditions, in atomic samples the effect 
of superradiant emission is observed^. A similar effect of 
coupling to common radiative modes appears also in the 
emission from QD systems. Here also superradiance-like 
phenomena occur in the evolution of the exciton occu- 
pation and polarization.— ~—. Apart from the modified 
evolution for spontaneous emission, collective coupling 
to the radiation field can lead to other physical effects, 



one of which is the vacuum induced coherence (VIC) 25 . 
This effect consists in spontaneous, partial coherent exci- 
tation transfer from an initially occupied QD to initially 
empty one which results in exciton occupation trapping 
in a decoherence-resistant state. Like the superradiance 
phenomena, the VIC effect can be expected to occur also 
in QD systems. 

On the other hand, although natural atoms and QDs 
have many properties in common, the characteristic fea- 
ture of the latter is the inhomogeneity of transition en- 
ergy which, for technologically realistic structures, is on 
the order of milli-electron- Volts. Moreover, small spatial 
separation of the dots leads to coupling between them^S. 
In our previous works we studied the impact of the tran- 
sition energy mismatch and the coupling between the 
two emitters on the stability of collective effect a 27 i 28 and 
their role in the linear— and nonlinear optical response 
of the DQD systems^. We showed that, in the absence 
of coupling the between the QDs, the appearance of an 
optically inactive subradiant state and a rapidly decay- 
ing superradiant one are extremely sensitive to the fun- 
damental transition energy mismatch and already for a 
mismatch on the order of micro-electron- Volts the col- 
lective character of the evolution is replaced by oscil- 
lation around the average exponential decay. The de- 
structive effect of the system inhomogeneity may be, to 
some extent, overcome by sufficiently strong coupling be- 
tween the dots which rebuilds collective behavior even in 
structures with technologically realistic values of energy 
mismatc h 22 ' 28 . We have also pointed out that phonon- 
induced dynamics can slow down the decay of a super- 
radiant state or speed up the emission from the super- 
radiant one^I. Based on these previous investigations of 
superradiance phenomena in QDs, one can expect that 
the VIC effect should, in principle, also be observable 
in these systems, at least in the presence of sufficiently 
strong coupling 27 . However, to our knowledge, its sta- 
bility against various inhomogeneities and perturbations 
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typical for the solid state environmnet (in particular non- 
parallel orientation of the interband dipoles in the two 
dots as well as phonon perturbations) have not been stud- 
ied. 

In this paper, we study the necessary conditions for 
the VIC effect to appear in a system of two vertically 
stacked semiconductor QDs. We investigate the role of 
the fundamental transition energy mismatch, coupling 
between the dots, and phonon-induced kinetics in this 
process. We also pay particular attention to the differ- 
ence of the magnitude of the interband dipole moments 
as well as to their non-parallel alignment (due to sub- 
band mixing). We show that in spite of all these inho- 
mogeneities that are typical to QD structures (as opposed 
to natural atoms) the vacuum-induced coherence can be 
almost fully stabilized in realistic pairs of non-identical 
QDs in a certain range of parameters. In particular, dif- 
ferent interband dipole moments for the two dots lead 
to the appearance of a state which is perfectly immune 
to radiative decay for a particular choice of the system 
parameters. 

The paper is organized as follows. In Sec. UH we 
describe the system under investigation and define its 
model. In Sec. IIII1 the method for describing the evolu- 
tion is described. Section |IV] contains the discussion of 
the results. Concluding remarks are contained in Section 
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II. THE SYSTEM 

The investigated system is composed of two vertically 
stacked semiconductor QDs interacting with quantum 
electromagnetic field and lattice vibrations. We restrict 
the discussion to the ground-level transitions with fixed 
polarization and spin orientations. We take into consid- 
eration only 'spatially direct' states in which electron- 
hole pairs reside in the same QD. Due to the strong 
Coulomb coupling these states have a much lower energy 
than the 'dissociated' states (the external electric fields 
which would change this pictur o 6 ' 32 are not considered in 
our discussion). In this manner, the DQD or QDM may 
be described as a four- level system, with the state 1 00) 
denoting empty dots, states 1 10) and 1 01) representing 
single-exciton states with electron-hole pairs residing in 
the lower or higher QD, and the 1 1 1) corresponding to 
the biexciton state, that is, to both QDs occupied with 
an exciton. 

As it is currently impossible to produce on demand 
pairs of QDs with identical fundamental transition ener- 
gies, we assume that the exciton transition energies for 
the two dots are different, 

Ei = E + A and E 2 = E - A, 

where E is the average transition energy and A is the 
energy mismatch. 

As in Ref. [29|, we describe the evolution in the 'rotating 



basis' defined by the unitary transformation 

U — e -i[B(|10)(10| + |01)(01|+2|ll)<ll|)+_H" rad +ff ph ]t/?i 

where -ff ra d and are the standard free photon and 
phonon Hamiltonians, respectively. The Hamiltonian of 
the system is then 

H = HuQD + #DQD-ph + -ffDQD-rad- (1) 

The first term describes exciton states in a DQD struc- 
ture. 

tfDQD = A(|10)(10| - |01)(01|) + Vb|ll)(H| 

+ V(|10}<01| + |01}<10|), (2) 

where Vb is a biexciton shift due to the interaction of 
static dipole moments and V is the amplitude of the cou- 
pling between the single-exciton states of the dots. 

The QDs are separated by a distance on the order of 
a few nm which is much smaller than the relevant wave 
length of the electromagnetic field with which the dots 
interact. This allows us to neglect the space dependence 
of the electromagnetic field and describe the coupling of 
excitons to the photon modes in the Dicke limit 20 . The 
relevant Hamiltonian in the dipole and rotating wave ap- 
proximations is then 

2 

^DQD-rad = £ £ ^X^ H, X + H.C 
a=l kX 

where <r" = are the creation and annihilation op- 

erators for the exciton in the ath QD, is the creation 
operator of the photon mode with the wave vector k, and 

V 2eQ€ r v 

is a coupling constant for the ath QD, where d a is the 
inter-band dipole moment for the ath QD, e\(k) is the 
unit polarization vector of the photon mode with polar- 
ization A, ojk is the corresponding frequency, eo is the 
vacuum dielectric constant, e r is the relative dielectric 
constant of the semiconductor and v is the normaliza- 
tion volume. We investigate only wide-gap semiconduc- 
tors with fundamental transition energies on the order of 
1 eV for which zero-temperature approximation for the 
electromagnetic modes may be used at any reasonable 
temperature. 

Interaction of the carriers confined in the two dots with 
phonon modes is modeled by the Hamiltonian 

^DQD-ph = £ C%&t £ f q a) {c q + d q ), 
a=l,2 q 

where c q and c q are creation and annihilation operators 

(1 2) 

of the phonon mode with the wave vector q and f q ' are 
the system reservoir coupling constants for the first and 
second QD, respectively. We model the electron and hole 
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wave functions by identical Gaussians with extensions I 
in the xy plane and l z along the growth direction, 



W(r) ~ exp 



1 x 2 + y 2 
2~ 



If! 

2 I 2 



For such wave functions and for the deformation potential 
couplings between the confined carriers and longitudinal 
phonon modes, the coupling constants have the form-l 



f(l,2) 
Jq 



± 



iq z D 



/ q exp 

where D is the distance between the dots and 



Hq 



2pvc\ 



exp 



Ijqj + l 2 q 



2„2 



Here u c / h are deformation potential constants for elec- 
trons/holes, p is the crystal density, c\ is the speed of 
longitudinal sound (linear phonon dispersion is assumed) 
and q± tZ are momentum components in the xy plane and 
along the z axis. 



III. THE EVOLUTION 

Analytical formulas describing the evolution of a pair 
of QDs are available for uncoupled systems (V = 0) in- 
teracting only with phonon modes 33 and, in the Markov 
limit, if only radiative decay is included 28 . In this paper 
we use the previously proposed method 3 ^ that allows us 
to deal with the simultaneous action of both the phonon 
and photon surroundings. Our approach is based on the 
equation of motion for the reduced density matrix of the 
carrier subsystem in the interaction picture, 



P = £rad[/>] + A>hM- 



(3) 



Here the first term accounts for the effects induced by the 
radiative environment, which is described in the Markov 
limit by the Lindblad dissipator 



£rad[p] = 
2 



afi=l 

where 



J2 I> ^(t)pai(t)-±{aP + (t)a«(t),p} 



(4) 



al(t)= [oS(t)] + =exp 



iH 



DQDt 



and 



Fa/3 — — 



E 3 



cr_ exp 



d a -de 



DQDt 



(5) 



If we assume that the spontaneous decay rates for the 
two QDs are Tn and r 2 2 then it follows directly from 
Eq. © that 
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Parameter 



Symbol Value 



Difference of deformation potential 

constants for electrons and holes 

Crystal density 

Speed of longitudinal sound 

Carrier localization extensions 

in the xy plane 

Carrier localization extensions 

in the growth direction 

Spatial separation of the dots 



G e — Oh 

Q 

Cl 



D 



9 meV 
5350 kg/m 3 
5150 m/s 

4.5 nm 

1 nm 
8 nm 



TABLE I. Parameters used in numerical simulations. The 
values correspond to a self- assembled InAs/GaAs system. 



where d a = d a /d a . By redefining the relative phase of 
the exciton states in the two dots one can assume without 
any loss of generality that Ti 2 and T 2 i are real. 

The second term in Eq. ([3]) accounts for the effects 
due to interaction with phonon surrounding, allowing for 
non-Markovian dynamics. To describe these effects we 
use the time-convolutionless equation 

£ P h[p] = 

dTlr ph [Hr,Qn-ph(t), [H D Q D ^ ph (r), p(t) ® p P h}]-, 

o 

where 



#DQD- P h(£) = exp 



iH 



DQDt 



#DQD- P h exp 



DQD£ 



p p h is the phonon density matrix at thermal equilibrium, 
and Tr p h denotes partial trace with respect to the phonon 
degrees of freedom. 



IV. RESULTS AND DISCUSSION 

Below, we present our results of simulations of the VIC 
process in pairs of QDs. In all the analyzed cases, we as- 
sume that the system is prepared initially in a localized 
state 1 10). In Sec. refsubseeddentical we explain the 
effect in a system of identical, uncoupled dots. Then, 
in Sec. IIVB1 we analyze the role of the relative magni- 
tude and orientation dipole moments and in Sec. IIV CI 
the effect of the energy mismatch and coupling between 
the dots. In Sec. IIVDI we analyze the interplay of all 
the parameters that distiguish QDs from natural atoms 
in the evolution of QDs interacting only with radiation 
reservoir. The phonon impact on the VIC is discussed in 
Sec. lIVEl 
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FIG. 1. (Color online) Vacuum-induced coherence in a system 
of two identical (A = 0, Tn = T22 = 1 ns -1 ) two-level sys- 
tems, (a) Exciton occupation of the system and occupations 
of the states 1 10) (nio) and 1 01) (noi). (b) The off-diagonal 
density matrix element poi.io- 



A. Identical QDs 

Collective coupling of two identical QDs (systems 
with identical fundamental transition energies and par- 
allel dipole moments of equal magnitudes) to the quan- 
tum electromagnetic vacuum leads to the appearance 
of a short-living (bright) superradiant state, |+) = 
( 1 10) + 1 01) ) /-\/2, and an optically inactive (dark) sub- 
radiant state |-) = (| 10) - 1 01)) /V2- The initial state 
of the analyzed system is a localized single-exciton state 
( 1 10) or 1 01) ) which can appear naturally, e.g., as an ef- 
fect of incoherent trapping or controlled tunnel injection 
of carriers in an injection structure similar to that stud- 
ied in Ref. HH. Such a state of a system of two identical 
and uncoupled QDs may be expressed as an equal combi- 
nation of the sub- and superradiant states, |— ) and |+), 



|10) = — (|+) + |-)) : 



|01> = 7! (| - 



-))■ (6) 



Coupling to the electromagnetic reservoir induces emis- 
sion from the bright state and, consequently, decay of a 
half of the initial excitation, while the other half is un- 
affected (trapped) due to the stability of the subradiant 
state for a system of two identical two-level systems^. 
This is observed as a spontaneous and coherent excita- 
tion transfer from initially occupied dot to the initially 
empty one until occupation of both QDs stabilizes at 
the same level (Fig. []Ji). While this process is taking 
place, coherence builds up spontaneously in the system 
(Fig.[Tp). As a consequence, the pair of identical QDs is 
trapped in a delocalized and decoherence-resistant state 
with the exciton occupation number equal 0.5 and real 
off-diagonal matrix element equal —0.25. This effect is 
referred to as VIC. 



B. The role of dipole moments 

Vertically stacked semiconductor QDs differ slightly 
in size and shape. If the system was formed in a self- 
assembled two-layer process then the upper QD is usu- 
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FIG. 2. (Color online) (a) Exciton occupation (red solid line) 
for a pair of decoupled (V = 0) QDs with identical transition 
energies (A = 0) but with non-identical values of the spon- 
taneous recombination rates (rn 7^ Tza)- (b-d) The total 
exciton occupation (red solid line) and the occupations of the 
individual dots (green and blue dashed lines) for the three 
cases shown in (a). The value of Tn = 1 ns -1 is the same for 
all the graphs in this figure. 



ally bigger than the lower one2&. The geometry and the 
structure of the QDs is reflected in the carrier wave func- 
tions and thus in the interband dipole moments. Accord- 
ing to Eq. ([5]), this leads to different spontaneous decay 
rates for the two QDs forming the investigated system. 
The dipoles corresponding to the upper and lower QD 
may differ in amplitude and, if the hole states in the two 
structures have different light-hole admixtures, also in 
orientation (see Appendix). 

The equation of motion describing the evolution of a 
system of two energetically identical ( A = 0) and uncou- 
pled (V — 0) QDs interacting only with the radiation 
reservoir is described by Eq. © with C p h — and £ ra d 
given by Eq. (Q} with cr±(t) = a±. It is known that 
in three-level open systems of this kind, non-radiating 
superposition states occur under certain conditions^! - — . 
In our case, a non-trivial (different than the ground state 
|00)(00|) stationary solution to the open system evolution 
equation, corresponding to the spontaneously formed, 
stable, delocalized dark state discussed above, exists only 
for Ti2 = VT1T2, i.e., for parallel dipole moments [Eq. 
(0] (such that di = d-i up to a phase factor). For dipoles 
that are parallel but have different amplitudes, the evo- 
lution is similar to the case of identical systems, i.e., cou- 
pling to photon surrounding leads to excitation transfer 
and occupation trapping. As can be seen in Fig. [2j the 
fraction of trapped exciton occupation depends on the 
values of the single dot decay rates and stabilizes at the 

level r 11 /(r 11 + r 22 ). 

As in the case of identical QDs, the suppression of 
the exciton decay is due to the existence of a dark state 
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FIG. 3. Exciton occupation of a pair of uncoupled (V = 0) 
QDs with identical transition energies (A = 0) and different 
spontaneous recombination rates (Fn 7^ T22) for a few values 
of the angle 9 between the dipole moments. The values of Tn 
and T22 shown in Fig. (a) are valid for both figures. 



FIG. 4. (a) The impact of the energy mismatch on exciton 
occupation of a pair of uncoupled QDs. (b) Exciton occu- 
pation of a technologically realistic QDM for a few values of 
coupling between the emitters (Fn = T22 = 1 ns _1 ). 



which, for T±i 7^ T22, is not strictly anti-symmetric, 



The role of the energy mismatch and coupling 
between the dots 



I dark) 



y/T n +T 



22 



(7) 



Due to unequal contribution from the single-exciton 
states to the dark state ([7]) the final occupation of the 
dots is non-symmetric, either, and equals 



and not 



rnr 2 



(T n + v 22 y 



respectively. Different dipole moments allow one to 
achieve different final situations. In the trapped state, 
the occupation of the initially excited dot may be lower 
[Figs.Etb) andCJc)] or higher [Fig.^d)] than that of the 
initially empty dot. 

Another factor that influences the interband dipole 
moments is the light-hole admixture. If this admixture 
is different for the two dots then the dipole moments 
corresponding to the QDs forming the QDM or DQD 
become non-parallel (see Appendix). This means that 
the off-diagonal decay rate r 12 = ^1^22 (l — # 2 /2) 7^ 
VT11I22, where 8 is the angle between the interband 
dipole moments of the two dots [see Eq. (IA.2I) ]. Conse- 
quently, the stationary solution to Eq. @ does not exist 
and quenching of the final exciton occupation is observed. 
The values of 9 are determined by the light hole admix- 
ture and typically are on the order of 0.01—. However, as 
can be seen in Fig. [3J even for much larger values of the 
angle between the dipoles, quenching of the occupation is 
weak. Although the system looses its coherence at long 
times, the character of the evolution remains the same as 
in the case of parallel dipoles on time scales much longer 
than the nominal exciton lifetime, i.e., initially excitation 
transfer between the dots takes place until the occupa- 
tion rate nio/noi close to that defined by the decay rates 
Tn and T 2 2 is reached and, then the impact of of hole 
subband mixing is manifested as a slow and equal decay 
of occupations [Fig. [3[b)]. 



In spite of rapid technological progress, manufacturing 
of DQDs or QDMs with identical fundamental transition 
energies is still not feasible. As can be seen in Fig.[4|a), 
for the realistic case of non-zero energy mismatch, the 
effect of trapping the system in an occupied and opti- 
cally inactive state is destroyed. Already for the energy 
splitting on the order of tens of /xeV quenching of the 
final occupation is observed (the relevant energy scale is 
the transition energy line width fiF = 1.7 /icV). QDs con- 
stituting such an inhomogeneous double dot systems in- 
teract with disjoint energy ranges of the electromagnetic 
field which destroys the collective character of coupling 
to photon modes. The energy mismatch of the dots slows 
the decay of the superradiant state and induces emission 
from the subradiant one2&. Due to the lack of a stable 
exciton state in which the system might be trapped the 
VIC effect in inhomogeneous QDMs is destroyed. 

Initially, the character of the evolution of an exciton oc- 
cupation of the inhomogeneous pairs of QDs does not dif- 
fer considerably from the corresponding case of a pair of 
identical dots. Until t ~ K/(2A), the coupling to photon 
modes maintains its collective character and excitation 
transfer from the initially occupied QD to the initially 
empty one takes place. Later, due to the emission from 
the subradiant as well as from the superradiant state, 
occupations of both dots decay^l. 

If the distance between the QDs is sufficiently small 
coupling between the systems (Forster or tunneling) be- 
comes effective and affects the evolution of carriers con- 
fined in the structure. Since sub- and superradiant states 
are eigenstates of the coupling part of the Hamiltonian 
[Eg. ((2J] separated by the energy 2V, sufficiently strong 
interaction between the dots rebuilds the collective char- 
acter of the interaction even in structures with techno- 
logically realistic energy mismatches on the order of 1 
me V 19 ' 28 . This also enables the VIC effect in inhomoge- 
neous DQDs to be rebuilt. As can be seen in Fig. 2Kb), for 
V A, a transition from a localized initial single-exciton 
state to a nearly stable state is possible for a system with 
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FIG. 5. Exciton occupation of a pair of technologically re- 
alistic QDs. (a) Total exciton occupation for A = 1 meV, 

V = —5 meV, Fn = 1 ns -1 , and for a few values of r 2 2- (b) 
The total occupation and the occupations of the individual 
dots for the last case shown in (a), (c) Exciton occupation 
for three sets of parameters stabilizing the occupation trap- 
ping. Red solid line: A = 0.5 meV, V = -4.47 meV, I\i = 
0.8 ns -1 , F 22 = 1 ns -1 . Green dashed line: A = 0.5 meV, 

V = -9.99 meV, Tn = 0.95 ns~\ r 22 = 1.05 ns 
dotted line: A = 1 meV, V = —5 meV, Tn = 1 ns 
r 22 == 1-49 ns -1 . (d) The total occupation and the occu- 
pations of the individual dots corresponding to the red line in 
(c). 



a technologically achievable energy mismatch. Although 
full stabilization of the VIC effect is impossible for non- 
identical dots and quenching of the exciton occupation 
always takes place sufficiently strong coupling between 
the dots considerably reduces the decay of the trapped 
state. For a weaker coupling between the dots, the exci- 
ton occupation decay is faster but still reduced compared 
to the system of uncoupled QDs. In contrast to identical 
dots, the two localized states have different contributions 
from the from sub- and superradiant states. This results 
in a different degree of trapping depending on the choice 
of the initially occupied stated. 



D. Interplay of energy mismatch, coupling between 
the dots and non-uniform dipoles 

For technologically realistic DQDs and QDMs, the fun- 
damental energy mismatch of the two QDs forming the 
system is on the order of milli-electron- Volts, the dipole 
moments differ slightly between the dots and the sys- 
tems are coupled with one another via Forster or tun- 
neling coupling. The evolution of such systems, for par- 
allel dipole moments, is shown in Fig. [5] In Fig. HJa), 
the dynamics of a system strongly stabilized by a cou- 
pling between the dots is presented for a few values of 
the decay rates. As can be seen, the rate of occupa- 
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FIG. 6. Exciton occupation of a pair of non-identical (A = 1 
meV, Tn r 22 ) and coupled (V = —5 meV) QDs for a few 
values of 8 (a) and for two values of F 22 (b). 



tion decay depends on the choice of the two single-dot 
decay rate parameters. In such a system, the coherent 
excitation transfer between the dots takes place on time 
scales similar to the previously discussed cases. However, 
the initial state is now a superposition of non-degenerate 
system eigenstates, which leads to fast oscillations of the 
exciton occupation between the two dots [Fig. |3£b)]. 

In the present general case, the spontaneously formed 
coherence and the trapped occupation can again be fully 
stabilized by an appropriate choice of parameters in a 
system with parallel dipoles. The single-exciton eigen- 
states of the system Hamiltonian (JTJ are 



|* 1 }=cos( I I |10)+sinl I 



101) 



(8a) 



1*5 



|10) 



|01), (8b) 



where tg(0) = V/A and tt/2 < <p < it/2. The corre- 
sponding eigenenergies will be denoted by Ei,E%. If the 
coupling between the single exciton states satisfies the 
relation 



V = - 



2Ayr77r 2 ^ 
|r 2 2-r u | 



(9) 



then one of the eigenstates [Eqs. (|8a|) and (|8b|) ] corre- 
sponds to the dark state ([7]). In such case, the evolution 
of a realistic system also leads to occupation trapping in a 
decoherence-resistant state [Fig. (He)] at the occupation 
I'll /(rn + T 2 2). As can be seen from Eq. ©, non- 
equal values of the spontaneous decay rate for the two 
dots open a possibility to stabilize the evolution of the 
system for wide range of parameters. As for a system of 
identical dots, the occupations of single dots also stabi- 
lizes, but again initial fast oscillations due to contribution 
from two energy eigenstates are observed [Fig. EJd)] . 

For the non-parallel dipole moments, i.e., for # 7^ 0, 
the effect of occupation trapping is in principle destroyed 
but for strongly coupled dots, the decay is very weak. For 
the interesting time scale on the order of ns the impact 
of non-parallel dipoles becomes visible only for the val- 
ues of 9 that exceed the realistic ones by an order of 
mag nitude^I (Fig.©. 
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FIG. 7. (Color online) The phonon spectral density for D — 
8 nm. 

E. The role of lattice dynamics 

Spontaneous emission from a system of QDs is af- 
fected by phonon dynamic o 31 i . Coupling between the 
QDs and lattice vibrations induces excitation transfer 
between the two single exciton eigenstates [Eq. (|5a|) and 
(|8b[) ]. For a pair of identical QDs (</> = ±7r/2), these 
eigenstates exactly coincide with the bright (superradi- 
ant) and dark (subradiant) states |+) and |— ) defined 
in Sec. IIV Al If the dots are non-identical, there is no 
perfect correspondence between these two pairs of states 
but, for non-zero coupling, still one of them is brighter 
and the other is darker. As we have shown previously*^, 
the phonon-induced redistribution of single exciton occu- 
pations among these two states strongly affects sponta- 
neous emission in a temperature-dependent way. For the 
system of vertically stacked dots analyzed in this paper, 
the amplitudes of both tunneling and Forster couplings 
are negative, so that cj> < 0. In this case, the darker 
eigenstate has a higher energy and will be affected by 
relaxation even at low temperatures. As we show below, 
phonon dynamics indeed breaks the relative stability of 
the darker state and leads to considerably accelerated 
decay of the excitonic occupation, except for special pa- 
rameter choices. 

While, in general, non-Markovian effects may be im- 
portant in the carrier-phonon dynamics that affects spon- 
taneous emission^, the phonon-related effects to be dis- 
cussed below can be understood within a Markovian pic- 
ture of transitions between energy eigenstates. The cor- 
responding rate for a transition between the eigenstates 
and |%) is 7i_>j = 27rsin 2 0J(|wy \)\n(wij) + 1|, 
where i,j = 1,2, toy = (Ej — Ei)/h and the spectral 
density is defined as 

9 

The spectral density for the inter-dot spacing D = 8 nm 
is plotted in Fig. [7] Apart from the usual cutoff at fre- 
quencies larger than l/c due to the restriction on the 
momentum non-conservation A(hq) < H/l, where I is the 
average dot size, it is characterized by the oscillations 
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FIG. 8. Phonon impact on vacuum-induced coherence in 
a pair of QDs with a fixed V/A ratio (corresponding to 
4> = —1.4): the comparison of the total exciton occupation 
in a system coupled to phonons and without phonons. (a) 
Q, — 4.0 meV; (b) short-time section of (a) with solid (red) line 
representing the total exciton occupation, dotted (magenta) 
line corresponding to the total occupation without phonons 
and the dashed (green) and dash-dotted (blue) lines showing 
the occupations of the two dots; (c) Q = 8.07 meV (corre- 
sponding to the 3rd minimum of the spectral density); (d) 
Q, — 13.36 meV (5th minimum of the spectral density). 



with the period Aw = 2nc/D which result from the dou- 
ble dot structure of the system^. 

The phonon-related effects on the spontaneous buildup 
of coherence is shown in Fig. [8j Here we keep the mixing 
angle tj> constant (there is, a constant ratio V/A) and use 
the energy splitting Ml = E\ — E 2 = 2y 'A 2 + V 2 as a 
parameter. In Fig.(5Ja) we plot the evolution of the total 
exciton occupation for Ml = 4.0 meV, where the spectral 
density is relatively large, corresponding to phonon tran- 
sition rates on the order of 1 ps _1 . One can see here 
that indeed the phonon-induced relaxation suppresses 
the vacuum-induced coherence, which is manifested by 
the rapid decay of the excitonic occupation. The role of 
phonons in this strong change of the occupation dynam- 
ics is clear in Fig. [5Jb), where we show the first 100 ps 
of the same evolution. A very fast phonon-induced redis- 
tribution of occupations of the two dots (shown in green 
dashed and blue dash-dotted lines) takes place on a pi- 
cosecond time scale after which the radiative effects are 
completely dominated by the phonon-induced thermal- 
ization which, at low temperatures, forces the system 
to stay in the bright state |\ta)- The detrimental ef- 
fect of carrier-phonon coupling can only be avoided if the 
phonon-induced occupation dynamics is made slow com- 
pared to the spontaneous emission. This can be achieved 
by using the oscillating form of the phonon spectral den- 
sity (see Fig. [7]). As can be seen in Fig.(SJc,d), If choosing 
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the energy splitting h£l such that it corresponds to the 
3rd minimum of the spectral density has little positive 
effect but for Ml in the 5th minimum, the phonon effects 
become very weak and the long living tail of the excitonic 
occupation is restored. 

An interesting additional effect can be seen if one com- 
pares Figs. Eta) and[H[d): If the phonon effects are strong 
then increasing the temperature slows down the occupa- 
tion decay. However, for weak phonon influence the tem- 
perature dependence is opposite. This can be explained 
by noting that in the presence of fast phonon-induced re- 
distribution of single exciton occupations, the occupation 
of these two states remains in quasi-equilibrium, which 
means that the occupation of the higher energy, darker 
state increases as the temperature grows. On the con- 
trary, if the phonon-induced dynamics is slow the sys- 
tem state is almost unperturbed and close to the spon- 
taneously formed dark superposition. Phonons lead to 
transitions out of this stable state with intensity growing 
with temperature. 



V. CONCLUSIONS 

We have studied the formation of vacuum-induced co- 
herence and the associated long-living trapped excitonic 
population in a pair of vertically stacked semiconductor 
QDs. We focused on the features that distinguish QD 
systems from natural atoms: the mismatch of transition 
energies, coupling, possibly non-identical dipole moments 
for the optical transition, and strong interactions with 
the phonon environment. 

We have shown that the VIC effect is very sensitive to 
the inhomogeneity of the QDs. Already for the funda- 
mental transition energy mismatch on the order of the 
emission line width, the exciton occupation is quenched. 
However, the destructive effect of the energy inhomo- 
geneity can be strongly reduced by coupling between the 
dots. While for dots with identical magnitudes of the 
interband dipole moments full stabilization of the spon- 
taneously formed coherence can only be achieved in the 
limit of infinite coupling, in pairs of QDs with different 
interband dipoles it is possible to adjust the energy mis- 
match, the coupling between the QDs, and the dipole mo- 
ments in such a way that a perfectly stable state forms 
in the spontaneous emission process and a fraction of 
the initial excitonic population attains a formally infi- 
nite life time. The non-parallel orientation of the inter- 
band dipoles, that may be caused by heavy-light hole 
mixing leads to negligibly weak effects in realistic struc- 
tures. While carrier-phonon coupling typically destroys 
the vacuum induced coherence on picosecond time scales, 
it can be overcome by appropriately selecting the energy 
splitting between the single exciton states. 

These results show that the VIC effect can be observed 
in realistic systems with energy splitting on the order of 
milli-electron- Volts provided that the system parameters 
(interband dipoles, coupling and energy mismatch) can 



be controlled with sufficient flexibility. This seems to 
be possible by appropriately designing the system on the 
manufacturing stage and then employing the dependence 
of various parameters on external fields (e.g, via Stark ef- 
fect or modification of electron-hole wave function over- 
lap). 

Let us note, finally, that the major experimentally 
detectable consequence of the appearance of vacuum- 
induced coherence in the double dot system is the long- 
living tail in the exciton occupation (or in luminescence 
intensity). Since this effect is of considerable amplitude 
and evolves on long, nanosecond time scales it should be 
relatively easily detectable with time-resolved lumines- 
cence or pump probe spectroscopy. 
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Appendix: Interband dipole moments and subband 
mixing 

In this appendix, we discuss the effect of light hole 
admixture on the non-parallel alignment of the interband 
dipoles for (nominally) heavy hole transitions. 

Within the 4-band Luttinger model^, the electron and 
hole wave functions are 

*L a W) = £u A (r, S )^(r) 

A 

*W(r )S )= M (r, S )^)(r), 

where a = 1,2 refers to a higher or lower QD, respec- 
tively, u\ is a Bloch function for the valence subband A 
(heavy holes: A = ±3/2, light holes: A = ±1/2), r is the 
space coordinate, s denotes spin, u c i is an electron Bloch 
function with a fixed spin s = 1/2 (for definiteness) and 
^e/h re f er to electron and hole envelope functions, respec- 
tively. 

The matrix element of the interband dipole moment 
for the ath QD is 

d a = Y,l d 3 r^\r, S )er^ c a \r,s). 

s J 

To calculate the above integral we sum over unit cells 
(labeled by R) and integrate over one unit cell (£ labels 
the position within the cell). The envelope functions vary 
slowly which allows us to assume that they are constant 
over one unit cell, ip e /h(R + £) « (p c /h(R). The Bloch 
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functions are periodic, s) — u(R, s), and orthog- 

onal for different bands. As a result, we can write^i 



A 

Here d x i is the bulk interband dipole moment, 



v is the unit cell volume, 

a[ a) = J d 3 R^(R)^\R) 



is the envelop function overlap integral, and we have re- 
placed the summation over unit cells with integration 
over R. 

We are investigating bright heavy hole excitons, hence 
a_ 3 / 2 ~ 1 and the other coefficients are much 

smaller. The non- vanishing bulk dipole moment matrix 
elements involving the spin- 1/2 electron state ar o 44 ' 45 



d_ 3 i = 

2 ■ 2 



— V3d* i i = — r=do 




Because the wave functions differ for the two QDs, the 
values of the heavy hole overlap integrals a^l may vary. 



Since, to the leading order, 
2 2 



Mo I 



» 



a_ 3 «o 



» 



» 



+ 



» 



» 
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(a) 

a_3 



the difference of the overlap integrals leads to different 
magnitudes of the dipole moments \d a \, a — 1,2. More- 
over, one has 



d\ ■ d% = a 



(1) (2)* 



Idol 



1 + 



(1) (2)* 
9 CI ifl ! 



„(1)J2)« 
1 a I a i 



Hence, 



d\ ■ d* 



\d x ■ d* 2 
Idillda 



1 



,(2) 



,(2) 



(2) 



(A.l) 



Thus, if the light hole admixture is different for the two 
dots, 



(1) (2) 

111 2 / 
(1) T 



a 



,(2) 



then tii • d\ 7^ 1, that is, the dipoles are non- parallel. 

As the subband mixing is typically small in self- 
assembled structures the angle 9 between the dipole mo- 
ments is small. Therefore, one can write 



di ■ d* 2 — e lv cost 



1-lt 



(A.2) 



where r\ is an irrelevant phase. Comparing Eqs. (|A.1|) 
and (|A.2|) one gets 
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